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ABSTRACT 

We present a cosmological hydrodynamical simulation of a representative volume 
of the Universe, as part of the Making Galaxies in a Cosmological Context (MaGICC) 
project. MaGICC uses a sub-grid feedback model that includes a thermal implemen- 
tation of both supernova -|- early (pre-supernova) stellar feedback that matches the 
observed stellar mass evolution of an galaxy at high resolution. The large volume 
simulation tests the stellar feedback at lower resolution across a wide range of galaxy 
masses, morphologies, environments and merging histories. 

The simulated sample compares well with a wide range of observed relations of 
high redshift galaxies (z ^ 2). It matches the slope and scatter in the stellar mass-halo 
mass (M^ — Mh) relation and the Galaxy Stellar Mass Function (GSMF) for low mass 
(M-t < 5 X lO^^M©) galaxies. The observed number density evolution of low mass 
galaxies is reproduced in our simulation, suggesting that our feedback implementation 
naturally detaches stellar mass growth from dark matter accretion history, which is a 
major improvement in galaxy formation models. The poor match of the M^, — M^^ and 
the GSMF at high masses (M^, ^ 5 x lO^^M© ) indicates that thermal feedback is not 
sufficient to limit star formation in these haloes. The simulations match the observed 
star formation history of the Universe and the star forming main sequence {SFR- 
M*). We see an increased specific star formation rate (sSFR) with redshift, in good 
agreement with observations that are corrected for nebular emission lines, although 
the simulations fall below the relation at z = 2. Altogether our results suggest that 
early stellar feedback, in conjunction with supernovae feedback, plays a major role in 
regulating the properties of low mass galaxies at high redshift. 

Key words: galaxies: formation - galaxies: evolution, interactions, structure - meth- 
ods: numerical ~ methods: N-body simulations 



1 INTRODUCTION 

Gravitational assembly of structure in a Lambda Cold Dark 
Matter (ACDM) Universe is well understood and mostly 
consistent with observations. However, the evolution of 
galaxies inside dark matter haloes presents many challenges 
for modellers. The baryonic physics in haloes is complicated. 
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involving various processes such as gas cooling, star forma- 
tion, radiative transfer, stellar and active galactic nucleus 
(AGN) feedback. These processes are highly non-linear and 
modelling them accurately is the major challenge for galaxy 
formation theory. 

Fortunately, there are now large catalogues of galac- 
tic data available at every epoch of cosmic history. They 
make it possible to compare galaxy formation models with 
observations, throughout their evolution. These catalogues 
include observations that present full spectral energy distri- 
butions of the galaxies from 1.4 Ghz radio continuum obser- 
vations with the VLA ( Karim et al.|20lT l, to infrared imag- 
ing with Hubble/WFC3, Spitzer/MIPS and VLT/HAWK-I 
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iKajisawa et al. 2010 Santini et al. 20121. These observa- 



tions give a complete picture of star formation even when it 
is dust obscured. 

From these observations, one can construct a cosmic 



star formation history to compare with models ( 


Lilly et al. 


1996'; 'Madau et al.|1996| |Hopkins|2004| |Wilkins ct al.|2008| 


Bouwens et al.||2012 


1. The shape of the cosmic star forma- 



tion history has a steep rise from z = to z = 1 before 
flattening off and then steadily decreasing from z — 2 to 
higher redshift. 

It is also possible to compare the star formation rate 
(SFR) of individual galaxies with their stellar mass (A/*) 
determined from infrared photometry. In observations, SFR 
and show a tight correlation that is sometimes called 



the star forming main sequence ( 


Brinchmann et al. 2004 


Noeske et al. 2007 Wuyts et al. 


20111. The slope of the 


relationship does not evolve significantly with redshift, but 



the normalization increases at higher redshifts (Whitaker 
et al.|2011||Kajisawa et al.|2010| l 



Dividing the SFR by M* gives the specific star forma- 
tion rate (sSFR), which provides a test of the star formation 
efficiency compared to prior star formation. Similar to the 
rise of the star forming main sequence, the sSFR rises with 



redshift (Karim et al. 2011 Kajisawa et al. 20101. Above 
z = 2, some observations show that the evolution of the 
sSFR flattens, though [Stark et al.| ( [2013[ ) found that when 
corrected for nebular line emission, the sSFR continues in- 
creasing up to 2; = 7. 

The primary constraint used for many models is the 
number density of galaxies as a function of their stellar 
mass, the galaxy stellar mass function (hereafter, GSMF). 
The GSMF is a Schecter type function characterized by a 
power law at low masses and an exponential cutoff. At z = 0, 



20091 



et al. (2012 



the exponential cutoff is at ~ 5 x 10 " Mq ( |Li fc White 
The GSMF evolves as a function of redshift: Santini 
I find that the low mass slope increases with red- 
shift while Peng et al. ( 20101) finds that the slope remains 
constant, but the normalization increases. 

Three types of models are commonly used to under- 
stand how stars populate galaxies: 

• Statistical models: compare statistics of simulations 
with observations 

• Semi-analytic models: populate dark matter haloes with 
stars based on halo mass, merger history, and single zone 
physics 

• Cosmological simulations: Model a volume of the Uni- 
verse with hydrodynamics 



1.1 Statistical Models 

The statistical models are based on comparing the GSMF 
with the dark matter halo mass function and lead to an 
understanding of how efficiently stars form as a function of 
dark matter halo mass. A set of cosmological parameters 



makes explicit predictions about the mass function (Press 



& Schechter 1974 Sheth et al. 2001 Reed et al. 2005 1 of 



dark matter haloes and how those haloes are distributed 
throughout the Universe. The observed GSMF has a differ- 
ent shape than the dark matter halo mass function found 
in simulations. The GSMF low mass slope (q) is shallower 
than the low mass dark matter mass function slope. The 



Mi, ~ 5 X 10^° M0 cutoff is at a lower mass than the dark 
matter mass function exponential cutoff. 

Halo Occupation Models make the reasonable assump- 
tion that the distribution of galaxies in the Universe is sim- 
ilar to the distribution of dark matter haloes, modulo some 
bias ( Peacock fc Smitli| |2000). Halo Occupation Models at- 
tempt to match the correlation function statistics of galax- 
ies and dark matter haloes to determine the stellar mass of 
galaxies that are most likely to be present in a particular 
dark matter halo. Using Halo Occupation Modelling, one 
can construct a Conditional Luminosity Function that can 



be compared with the real luminosity function (Yang et al. 
van den Bosch et al.|2007"| ) 



2003 



Conroy et al. ( 2006 1 realized that if one used satel- 
lite masses at their time of accretion, then the clustering 
statistics of mass-ordered dark matter halo samples matches 
galaxies. This realization lead to the abundance matching 
technique in which galaxies are placed in dark matter haloes 
with the same stellar mass ranking as that of the dark mat- 



ter halo mass rank ( Conroy fc Wechsler|2009 Moster et al. 
20T0| |Guo et al.||2010| [Behroozi et al.poTop Such a match 
leads to the stellar mass-halo mass (M« — Mh) relationship, 
the key constraint for our model. The — Mh relation con- 
sists of two power laws with a steep slope at low masses and 
shallow slope at high masses. Stars form most efhciently at 
the break mass. The star formation efficiency drops quickly 
to both higher and lower masses. The characteristic mass of 
the break in the power law is Mhaio ~ lO^^M© at 2: = 
( Moster et al]|2 013). 

The availability of luminosity functions at high redshifts 
means that we can trace the evolution of the — Mh rela- 
tion. Abundance matching indicates that the M« — Mh re- 
lation evolves surprisingly little ( [Behroozi et al.|[2013[ ). The 
star formation efficiency evolves most significantly to higher 
halo masses at higher redshift from Mh ~ 10^^ Mq at z = 
to Mh ~ 1 0^^ M0 at z = 3 ( [Moster et "aLl[2013|Behroozi 
et al.|2012| ). 

The key finding of the abundance matching models is 
that the star formation peaks earliest in the highest mass 
galaxies whereas, in the lowest mass galaxies, the SFR in- 
creases monotonically with time. This is a reflection of 
galaxy downsizing. This represents a delay of star formation 
in low mass halos and is the most important feature that 
must be reproduced in models in order to get the evolution 
of low mass galaxies right. 

1.2 Semi- Analytic Models 

Semi-Analytic Models (SAMs) try to match the GSMF at 
z = G using physical prescriptions based on the mass and 
merger history of dark matter haloes taken from simula- 
tions ( [Kauffmann et al.|1993[ ). SAMs show that supernovae 
can limit star formation in low mass galaxies ( White & Rees 



T978[|White fc Frenk|1991||Somerville fc Primack|1999||Ben 
son et al.[|2003| ) and that active galactic nuclei (AGN) can 
limit star formation in high mass galaxies (e.g. De Lucia 
et al.r2006. Bo wer et aL[[2006[ ) 



While SAMs do well matching the evolution of the high 
mass luminosity function, they do not match the evolution 
of low mass galaxies at high redshift ( Guo et al.|2011 1 . Cur- 
rent SAMs include strong stellar feedback to reproduce the 



GSMF at z = (e.g. [Guo et al.[[2011[ [Bower et aL][2012[ ) 
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but the low and intermediate mass galaxies build their stel- 
lar mass at early times {z > 2) following the assembly of 
the dark matter mass, because the feedback mechnism in 
these SAMs do not delay star formation in low mass halos. 
That means that there is little evolution in the SAM lumi- 
nosity functions after 2; = 2, in contrast with observations 
(e.g. Fontanot et al.|2009||Marchesini et al.,2009,|Guo et aL] 



and many other observed galaxy properties at « = ( Oppen 



2011 



The early star formation means the SAM galaxies 



have low specific star formation rates at z < 2 (e.g. Daddi 



et al. 2007 Damen et al. 2009 1 and high values at z > 3 
(e.g. Bouche et al.|2010|Dutton et al.|20To Weinmann et al 



2011) although this picture might change at high redshifts 
due to refinement in the observational estimates of sSFR 
I Stark et al.||2013 L This discrepancy has been looked at in 
detail by Weinmann et al.] p012), who use the number den- 
sity evolution of low mass (9.27 < log{AU/MQ) < 9.77) 
galaxies as a diagnostic to find that the observed evolution 
of the number density is not reproduced in any SAMs or 
simulations. They argue that the simple supernova feedback 
mechanism used in these models that gets the present day 
GSMF correct does not decouple star formation from the 
parent DM halo growth. 

1.3 Simulations 

Hydrodynamical simulations model galaxy formation with 
more resolution than SAMs and include self-consistent in- 
teraction of dark matter and baryon evolution. Although 
the efficiency of computational calculations has increased, 
it is still not possible to resolve many important physical 
processes, so they must be modelled at the 'sub-grid' level. 
These processes include gas cooling, star formation and stel- 
lar feedback. 

Since relatively little is known about star formation 
and feedback, the models include free parameters, which 
are constrained by observations. Star formation model pa- 
rameters are constrained using local observations of the 
Kennicut-Schmidt gas density-star formation density rela- 
tion ( Springel fc Hernquist|2003 Stinson et al.|2006 Schaye 
|fc Dalla Vecchia||2008| ). The energy feedback from stars is 
modelled either by adding velocity to gas, called kinetic feed- 
back, or adding thermal energy and delaying cooling, called 
thermal feedback. These models have been constrained based 



on observations ( 


Springel & Hernquist 


2003 


Oppenheimer 


& Dave"2006| |Dalla Vecchia & Schaye||2008 


|Crain et al. 


2009; McCarthy et al.||2012|. The model used 


in this paper 



instead constrains stellar feedback to match the evolution of 
Mi, - Mh relation. 

There has been a lot of research on the optimal velocity 
for winds driven using kinetic feedback. The original mod- 



els used a fixed wind velocity (Springel & Hernquist 2003 
[Grain et al.|2009| [McCarthy et al.|2012| |, however, they had 
difficulties reproducing the GSMF at z = 0. Observations 
of metal absorption lines in outflows show that wind ve- 
locities are not constant, but are correl ated with star for- 
mation rate, « ^ ^ q ([Martin[[2005[) and 



SFR"-' at z ^ 1.4 (Weiner et al. 20091. These obser- 



vations motivated using momentum conserving wind mod- 
els in which mass loading depends o n the mass of the host 



galaxy such that M^ind/Mi, oc V^^^^ (Oppenheimer & Dave 
2006} [Oppenheimer fc Dave|2008[ [Dave et al.|2011a|b[ ). Mo- 
mentum conserving winds successfully reproduce the GSMF 



heimer et al.|2010|Dave et al.|2011a|b[;[Puchwein fc Springe! 



20131, but has similar shortcomings with the low mass end 



of the luminosity function as the SAMs at high redshift. 



Weinmann et al. (20121 conclude that the current models 
of stellar feedback (in both SAMs and simulations) are un- 
likely to decouple the galaxy and DM halo growth due to its 
fundamental dependence on host halo mass and accretion 
history. 

In thermal feedback, stars heat the surrounding gas par- 
ticles adiabatically, which builds up pressure that can push 



gas out of galaxies ( Gerritsen fc Icke 1997 Thacker & Gouch- 



2000 Kawata fc Gibson 2003 Stinson et al. 20061 



Since stars form in very dense regions, the cooling time of 
the surrounding gas is short, making it necessary to account 
for unresolved inhomogeneities in the inter stellar medium. 
One method to account for such 'sub-grid' physics is to delay 
cooling within the blast radius of supernovae. Simulations 
using thermal feedback have so far focused on disk structure 



2010 



[20TT 



using high resolution zoom in simulations ( Governato et al. 
[Brook et aLlpOTT] [Sawala et aL][2011[ [Guedes et al. 
Agertz et aLpOlTf ^ Some recent simulations of a hand- 



ful of galaxies have indicated that adiabatic feedback pro- 
duces galaxies that follow — below Mhaio < W^^Mq 
(Brook ct al.|2012[ [Mimshi et al.|20T2 |. 

In most models of stellar feedback, only feedback from 



supernovae is considered, but Murray et al. (20101 recog- 



nized the amount stars can disrupt molecular clouds be- 



fore any stars explode as supernovae. Hopkins et al. (20111 



and Agertz et al. (20121 implemented early stellar feedback 
schemes that rely on IR radiation pressure and tested them 



on isolated galaxy simulations. Lopez et al. (20111 and Pel 



legrini et al. (2011 1 found that when they mapped the pres- 
sure in different phases of the gas in the 30 Doradus region 
of the LMG, UV photoheating provides more pressure than 
IR radiation pressure. 



In Stinson et al. (20131, we assume that photo-heating 
from massive stars is thermalised by the time it reaches the 
spatial scales resolved in cosmological simulations. So, we in- 
ject thermal energy . equal to the fraction of the bolometric 
luminosity emitted in the UV in the time between the forma- 
tion of the star and the first supernova explosion. This early 
stellar feedback limits star formation to the amount pre- 
scribed by the M* — Mh relationship and delays star forma- 
tion in an L* galaxy, so that the galaxy follows the evolution 
of the - Mh relationship. This is a major improvement 
over previous galaxy formation models, as the delayed star 
formation means that star formation is decoupled from DM 
halo mass growth. Some side-effects of using early stellar 
feedback include transforming DM cusps to cores in galax- 
ies up to L* masses (Maccio et al. 20121 and populating 



the circum-galactic medium with hot metal enriched gas, 
matching OVI observations ( Stinson et al.[[2"013 l. 

In this paper, we explore how the early stellar feedback 
model, described in Stinson et al. (20131, affects the global 



properties of galaxies on a large scale. To study this we sim- 
ulate a large volume of the Universe, 114 Mpc on a side, 
as part of the Making Galaxies in a Gosmological Context 
(MaGICC) project. This simulation tests the effectiveness 
of our model at low resolution across a wide range of galaxy 
masses, environments and merger histories. We compare the 
properties of the galaxies in our simulations with observed 
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statistical properties of high redshift galaxies like the GSMF, 
stellar to halo mass relationship, star formation rate, and the 
number density evolution of low mass galaxies through cos- 
mic time. In Section [2] we briefly outline the star formation 
and stellar feedback mechanisms used in our simulations, 
in Section |3] we present our results at z ^ 2 and compare 
them to the current observational estimates. In Section|4]we 
summarize our results and discuss future challenges. 



2 SIMULATION METHOD 

We simulate a cosmological volume, 114 Mpc on a side, from 
2 = 99 to z = 2. It is created using WMAP7 initial condi- 
tions with (ft, Qm , ^A, Ob , 0-8 ) = (0.702, 0.2748, 7252, 
0.0458, 0.816) ( [Larson et al.|2011|p<omatsu et al.|2011| . The 
simulation includes 512^ dark matter and 512'^ gas particles. 
The dark matter particles each have masses of 3.4 x 10* 
M0and a softening length of ~ 3.7 kpc. Initially, gas par- 
ticles each have masses of 6.9 x 10^ Moand star particles 
initially have masses of 1.3 x 10^ M©. Gas and star particles 
have a softening length of ~ 2.17 kpc. 

All the simulations use the smoothed particle hydrody- 
namics (SPH) code GASOLINE ( Wadsley et al.|2004 l. Details 
of the physics used in the MaGICC project are detailed in 



Stinson et al. ( 2013 1. Briefly, stars are formed from gas cooler 
than Tmax = 10 K, and denser than 9.6 cm~^ according to 
the Kennicutt Schmidt Law as described in IStinson et al.l 
(20061) with c*=0.1. 



The star particles are massive enough to represent an 
entire stellar population consisting of stars with masses 
given by the Chabrier (20031 initial mass function. 20% 
of these have masses greater than 8 Mq and explode as 
Type II supernovae from 4 until 35 Myr after the stellar 
population forms according to the Padova stellar lifetimes 
( |Alongi et a"L]p93l [Bressan et al] [19931 ). Each supernova 
ejects EsN = ergs of purely thermal energy into the 
surrounding gas (~ 1 kpc at the resolution of our simula- 
tions). The supernova energy would be radiated away before 
it had any dynamical impact because of the high density of 
the star forming gas. Thus, the supernova feedback relies on 
delaying the cooling based on the sub-grid approximation of 



a blast wave as described in Stinson et al. ( 2006 1 



The supernovae feedback does not start until 3.5 Myr 
after the first massive star forms. However, nearby molec- 
ular clouds show evidence of being blown apart before any 
SNI I exploded ([Murray et al.|20lT| and [Lopez et aL\ ( [2011 1 
and Pellegrini et al. (2011) found that most of the pres- 



sure comes from UV photoheating is the dominant feedback 
mechanism in early phases of star formation by mapping out 
the pressure in different phases of the gas. In simulations in 
the MaGICC project, like those here, 10% the UV lumi- 
nosity of the stars is injected into the surrounding gas over 



this 3.5 Myr period without disabling the cooling. (Stinson 



et al. 2013 1 showed that this energy limits star formation 
to the amount prescribed by the M+ — Mh relationship at 
all redshifts. The current work is our attempt explore how 
this star formation and feedback prescription works at lower 
resolutions over a wide range of galaxy masses. 



2.1 Halo identification 

For each snapshot, we find all the virialized haloes within the 
high resolution region using a Spherical Overdensity (SO) 
algorithm. Candidate groups with a minimum of Nf = 100 
particles are selected using a FoF algorithm with linking 
length (j> = 0.2d ~ 22 kpc (d is the average particle sepa- 
ration). We then: (i) find the point C where the gravita- 
tional potential is a minimum; (ii) determine the radius f of 
a sphere cantered on C, where the density contrast is Avir, 
with respect to the critical density of the Universe. Using all 
particles in the corresponding sphere, we iterate the above 
procedure until we converge onto a stable particle set. We 
use a time varying virial density contrast Avir, whose value 
is determined using the fitting formula presented in |Mainiiii[ 



et al. (20031, based on linear theory and the spherical top- 



hat collapse model. We include in the halo catalogue all the 
haloes with more than 100 particles (see Maccio et al. [20(37 



[2008| for further details on our halo finding procedure 



3 RESULTS 

We compare the simulated galaxy population in a 114^ 
Mpc'^ volume with a set of basic properties derived from 
the most recent observational estimates. These include the 
galaxy stellar mass function (GSMF), stellar mass-halo 
mass (Af* — Mh) relationship, cosmic star formation history 
(SFH) , star forming main sequence, and specific star forma- 
tion rates (sSFRs). Individual galaxies have been shown to 
match observations well (Brook et al. 2012 Stinson et al 
I2013|, so the volume provides an opportunity to test the ac- 
curacy and effectiveness of this feedback model at low reso- 
lution and high redshift. All the observational estimates of 
stellar masses and SFRs, that our results have been matched 



to, have been corrected to a Chabrier ( 2003 \ IMF. 



3.1 Stellar - halo mass (M* - Mh) relation 

Figure [l] shows the M* — Mh relation for all the galaxies 
in the simulated volume (black points) that contain a mini- 
mum of 20 star particles, or a stellar mass of ~ 3 x IO^Mq. 
The galaxies trace the slope of the M* — Mh (green line) up 
to Mhaio = IO^'^Mq at all redshifts where it has been ex- 
amined. The scatter of the simulated galaxies also matches 
the variation in the relation as obtained by [Moster et al.| 
(20131 (grey shaded area). The agreement points to the fact 
that the stellar feedback effectively regulates star formation 
to produce the right amount of stellar mass in a given halo 
mass at all times. 

Above a halo mass of 10^^ M0, abundance matching 
(green line) shows a decrease in star formation efficiency. 
This is not reproduced in the simulation. The star forma- 
tion efficiency actually increases at Mhaio ~ 4 x IO^'^Mq, 
before decreasing slightly as represented by the slightly shal- 
lower slope of the simulation points. The reduced SFE is 
due to the reduced gas accretion because of the high virial 
gas temperature of the halo. However, this slight decrease 
in SFE does not reduce the star formation in these high 
mass haloes to the extent observed. The implemented stel- 
lar feedback model is insufficient in these high mass haloes. 
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Figure 1. M^, — relation at different redshifts. The blue points are simulated galaxies, the green line is the [Moster et aL] | |2013| l 
relation derived from abundance matching techniques and the grey shaded area is the scatter derived for the relation. Our simulated 
galaxies match the relation below M/^aio < IO^^Mq, but star formation is too efficient in high mass haloes. 



Some other quenching mechanism is required such as feed- 
back from a central super-massive black hole (AGN feed- 
back, e.g. Fanidakis et al.pOll Springel et al.|2005l. 



whereas the observations show an exponential cutoff. This 
again indicates that stellar feedback is insufficient to limit 
star formation in these high mass haloes. 



3.2 The galaxy stellar mass function (GSMF) 

The GSMF measures the number of galaxies of a certain 
stellar mass in a given volume of the Universe. The era of 
deep, high redshift surveys has provided detailed GSMFs 
out to 2 = 3. We compare our simulation results to ISantim] 



et al. (20121, who use deep WFC3 near-IR data and deep 



Hawk-I Ks band data, to determine the low-mass end of 
the GSMF. The observed GSMFs are presented for various 
redshift ranges. To compare with them, we use the simu- 
lated GSMF from the middle of the observed redshift range. 
Figure [3] shows that the simulated galaxies (blue line) trace 
the low mass end of the observed GSMF (red points) very 
well. The feedback model makes the slope of the GSMF 
as shallow as the observed value, which is non-trivial and 
is a major improvement over previous attempts to match 
the GSMF at high redshift (e.g. |Guo et~aL][25dlt . A smaU 
discrepancy remains, as the simulated number density of 
high mass galaxies continues to decrease at the same rate. 



3.3 Number density evolution of low mass 
galaxies 



Weinmann et al. ( 2012 1 used the number density evolution of 
low mass (9.27 < log(M*/M0) < 9.77) galaxies to show that 
semi-analytic models or cosmological hydrodynamic simula- 
tions do not correctly model low mass galaxies. They argue 
that the simple supernova feedback mechanism changes the 
stellar mass at z = 0, but renormalizes the star formation 
history and thus does not decouple star formation from DM 
accretion. Stinson et al. ( 2013 1 showed for a single high reso- 



lution L* simulation that early stellar feedback can break the 
coupling of star formation to dark matter accretion. Figure 
[3] shows the number density evolution of low mass galaxies 
at high redshifts in our simulation volume compared to ob- 



servations taken from Figure 1 of Weinmann et al. ( 2012 1 , as 



weU with the SAM described in Guo et al. ( 2011 1. The simu- 



lation matches the observational results much better and lies 
well below the values obtained by the SAM. We note that 
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9.0 9.5 10.0 10.5 11.0 11.5 12.0 9.0 9.5 10.0 10.5 11.0 11.5 12.0 



Figure 2. Galaxy stellar mass function atz = 2& 2 = 3 (blue lines). The observational data points (rod points) are taken from |Santini| 
|et al.| | |2012l l. The low mass slope of the GSMF matches with observations, but there is no exponential decrease in the high mass end. 



Total SFR 

5ra-9.5 <log{M,jM^) < -10.5 

SFR -10.5<(fjff(Jtf,/jWg)<=11.5 




Figure 3. Number density of low mass galaxies as a function of 
redshift (black line) compared to various observational estimates 
(points) and a semi analytic model (green line) as described in 
|Guo et al.| | |20lT| l. Our simulation results go through the obser- 
vational points, but the slope is steeper. 

the slope obtained from our model is still slightly steeper 
than observed, indicating that the simulation is building low 
mass galaxies faster than observed. 



Figure 4. The evolution of the star formation rate density. The 
blue points are a compilation of star formation rate density es- 
timates taken from |Moster et al.| ( |20l"3] l. The black solid line is 
our result for all galaxies, the blue line is the SFR density for low 
mass galaxies (9.5 < log(M*/M0) < 10.5), the green for inter- 
mediate (10.5 < log(M*/M0) < 11.5) mass galaxies and the red 
line for high mass (log(M*/M0) > 11.5) galaxies. There is a clear 
trend of decreasing star formation at ^ ^ 3.5 in the highest mass 
galaxies. 



3.4 Star formation History 

We can also compare our simulation with the total number 
of stars formed in the Universe as a function of time. Fig- 
ure ID shows how the cosmic star formation rate evolves as 
a function of redshift ('Lilly-Madau plot') in our simulated 
volume. The observed points used for comparison are taken 



from Moster et al. (20131 and include star formation esti- 



mates derived from rest frame UV (|Salim et al.|2007| van der 


Burg et al.||2010||Robotham & Driver||2011| Bouwens et al. 


2011 Cucciati et al.||2012l, Hq (Ly et al.||2011 


1, combined 


UV and IR ( Zheng et al. |2007 Kajisawa et al. 


20101, FIR 


(Rujopakarn et al. 2010 


1 and radio observations (Smolcic 


et al.||2009 Dunne et al. 


2009 Karim et al.||2011 1. The to- 
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tal SFR density (black line) passes through the observations 
from 2 = 2 — 5. 

The total SFH can be divided into separate lines based 
on the stellar mass of the halo at 2; = 2 in which the stars 
are formed. The lowest mass galaxies (9.5 < log(M*/M0) < 
10.5) contribute little to the overall SFR density, while the 
the intermediate (10.5 < log(M*/MQ) < 11.5) and high 
mass (log(M*/MQ) > 11.5) contribute equally up to z = 
3. Below this redshift, the SFR flattens out in the highest 
mass galaxies. This flattening is not sufflcient to explain the 
quenching of high mass galaxies as shown by the failure of 
the simulated A/* — Afh and GSMF relations at the high mass 
end. We note that our match of the star formation history is 
not greatly affected by the excess star formation in galaxies 
with (log(M*/M0) > 11.5) because even though the galaxies 
in that mass range form too many stars at z 5C 3.5, they are 
not the dominant population of galaxies at those redshifts. 



Simulations 

i stark etal. 2012 



7 
en 

Ll. 

to 

lo" 



3.5 Star forming main sequence 

Observations show that star-forming galaxies have a tight 



correlation between their SFR and M* (e.g., 


Elbaz et al. 


2007 


IPannella et al.|2009||Wuyts et al.|2011||W 


litaker et al. 


2012 


I. This correlation has been called the "star forming 



mam sequence. 

We compare the SFRs of our simulated galaxies with ob- 



servational estimates by Kajisawa et al. ( 2010 1 and Whitaker 



et al.| ( |2012[ ). |Kajisawa et al.| ( |2010| ) studied SFR as a func- 
tion of M« for galaxies at 0.5 < z < 3.5 in the GOODS- 
North field, using the K-selected sample from Subaru- 
MOIRCS. They determined SFRs from rest-frame, dust- 
corrected UV luminosity and the Spitzer-MlPS 24 fim flux. 
The depth of their data allowed them to constrain the slope 
of the SFR-M* relation down to M^, = lO'^'^M© at z ~ 3. 
The median SFR as a function of stellar mass (green curve) 
from their sample of galaxies is plotted in top panels Fig. [5] 
at z = 2 & 3. The slope of their relation is close to unity 
for low mass galaxies at these high redshifts. Our simulated 
galaxies match these observations well at z = 3, but have 
nearly two times less star formation at 2: = 2. This dis- 
crepancy at 2 = 2 presents a challenge for all hydrodynamic 
simulations and SAMs ( [Weinmann et al.|20lT| ). |Dave|f2008| ) 
suggested that an evolving stellar IMF is required to reduce 
the discrepancy in this relation out to z = 2. 

Whitaker ct al. (201^ measure star formation rates us- 
ing the NEWFIRM Medium-Band Survey from MIPS 24 fim 
fluxes. At 2: > 2 their detection limit is Iog(M*/M0) > 10.7. 
For these galaxies, they find a shallower, sub-linear, slope 
for their star forming main sequence, SFR oc M" "**, with 
a constant scatter of 0.34 dex. Above their detection limit, 
our simulated galaxies (black points) match the observations 
(red line) well as seen in the top left panel of Fig. [5] 



3.6 Specific star formation rate evolution 

Another common way to compare star formation rates with 
galaxy stellar masses is the specific star formation rate 
(sSFR), which gives the amount of star formation in haloes 
for a unit stellar mass of material. As one would expect 
from the star forming main sequence, the bottom panels of 
Fig.jsjshow that the simulated galaxies (black points) match 



Figure 6. The evolution of the specific star formation rate for 
sample galaxies which have stellar masses within a narrow range 
around ~ 5 X 10^ Mq (blackline). We compare the simulations 
with the observational estimates of [Stark et al.| l |2013| . The evo- 
lution of sSFR matches well at z ^ 3, but is below the observed 
value at lower redshifts. 



the Kajisawa et al. (20101 (green curve) observed sSFRs at 
z = 3 but have ~ 2 times lower sSFR at z = 2. We also 
compare our simulated results with 1.4 GHz radio contin- 



uum observations from Karim et al. (20111 of star forma- 
tion in galaxies in the 2 deg COSMOS field. The simulated 
galaxies in our volume (black points) are in good agreement 
above log(Mj,/MQ) > 10.7, but are 2 — 3 times lower below 
this mass range at both z = 2 & 3. 



Karim et al. 


(2011), like other authors before them 


( Stark et al. 2009 


Gonzalez et al. 2010 1, found that sSFR in- 



creases for galaxies in a given stellar mass range from z = 
to z ~ 2, but then does not evolve much from z = 2 to 



7. Weinmann et al. (20111 shows that such observa- 



tions are contradictory with most models in which higher 
gas accretion rates at higher redshift and lower galaxy stel- 
lar masses translate into larger sSFR in galaxies within a 
fixed stellar mass range. 



Stark et al. (20131 re-examined their data and found 
that their Spitzer-IRAC photometry was contaminated by 
nebular emission. They use the photometric excesses in the 
contaminated [3.6] filter to estimate the equivalent width 
distribution of Ha emission at 3.8 < z < 5.0. The corrected 
sSFRs increase from z = 4toz = 7bya factor of ~ 5 
similar to what happens model predictions. Figure |6] shows 
the evolution of the sSFR in simulated galaxies (black line) 
within a narrow stellar mass range around ~ 5 x lO^M© . The 
simulation values are consistent with the corrected IStarkI 
et al. (20131 values (green points) for z > 3. However at 
z < 3 our simulation results are below the observed relation. 

Many other authors find also find lower than observed 
sSFRs in their models ( |Dave||2008| [Weinmann et al.||2012 |. 
The higher observed sSFRs again indicates delayed star for- 
mation in low mass galaxies. Although our model does a 
better job of delaying the star formation at early times than 
most SAMs and hydrodynamic simulations, below z — 3 
the simulated haloes may still be forming too few stars. 
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Figure 5. Top panels: The star forming main sequence (black points) at different redshifts. The 2 = 2 result for galaxies with 
log(M*/MQ) > 10.7 is matched to the observational results of Whitaker et al.] | |2012| (Red line). The slope of the main sequence 
is much steeper at lower masses. This matches well with the observational estimates derived by |Kajisawa et al.|||2010|| at 2 = 2fc3 ( green 
curve). Bottom panels: The simulated sSFR (black points) matched to observational results matched from |Karim et al.| l [2011| (Red 
points) anc |Kajisawa et al.| (j20ToJ (green curve) . The simulated sSFR lies below the observed values for low mass galaxies at z = 2 by a 
factor of ~ 2, but matches very well at z = 3. 



This suggests the importance of some other physical mecha- 
nism, not modelled in our simulation, like the dependence of 
the star formation on gas metallicity ( Krumholz et al.|2011 



Krumholz fc Dekel||2012 1, that could further delay star for- 
mation at earlier times and increase the sSFR of these galax- 
ies at z = 2 



4 DISCUSSION AND CONCLUSIONS 

We examine the effect of early stellar feedback used in 
the Making Galaxies in a Cosmological Context (MaGICC) 
project on a broad sample of galaxies in a cosmological vol- 
ume of 114"^ Mpc"^. The stellar feedback used is exactly the 
same as that used for a high resolution galaxy. We com- 
pare the simulated galaxies with the observed — Mh rela- 
tion, the galaxy stellar mass function, the cosmic star forma- 
tion history, the star forming main sequence and the specific 
star formation rate. The simulated galaxies do a good job 
matching each observation to z = 2, the time when previous 
models have most deviated from observations. Our use of 
thermal stellar feedback is the key difference between our 
simulation and ones run previously. The way that it delays 
star formation in Mh < IO^^Mq galaxies allows the simula- 
tions to match the observations. 

Our simulation does not yet follow the evolution of the 
galaxy population all the way to z = 0, but we think the 
results at 2 = 2 are significant enough to report and it will 
take a significant amount of time for the simulation to reach 



z = 0. We are encouraged by our previous simulations of 
the evolution of individual high resolution galaxies using the 
same stellar feedback match the evolution of an L* galaxy 
( Brook et al.|2012| [sFinson et al.|2013| |. It seems hkely that 
our simulated low resolution galaxies will also continue to 
match the low mass end of the galaxy stellar mass func- 
tion at z = 0. We plan to address this point in detail in a 
forthcoming paper. 

The simulated galaxies not only follow the Af* — Alh for 
Mh < IO^^Mq at all the redshifts examined but also match 
the scatter in the relation. Correspondingly, the simulated 
galaxies match the shallow slope at the low mass end of the 
galaxy stellar mass function. The slope of the GSMF rela- 
tionship was not a constraint for the simulation, but is a 
natural by product of the stellar feedback recipe used. It is 
a major improvement over previous attempts to match the 
GSMF at high redshift. The early stellar feedback decouples 
the growth of stellar mass from DM mass by effectively blow- 
ing the gas away from the disc either into the circum-galactic 
medium or entirely out of the halo. This helps regulate the 
number density of low mass galaxies to the observed values 
by delaying star formation in these haloes. 

The simulated star formation history of the Universe 
also matches a variety of different observations. The model 
predicts that the lowest mass galaxies (9.5 < log(M*/Mo) < 
10.5) contribute little to the overall SFR density, while the 
intermediate (10.5 < log(A/*/M0) < 11.5) and high mass 
(log(Af*/M0) > 11.5) galaxies contribute equally up to 2 = 
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3. After z — 3, the star formation slows in the highest mass 
galaxies. 

At Alh > IO'^^Mq, too many stars form, which is shown 
by the presence of galaxies above the abundance matching 
M* — Mh relation and the lack of an exponential cutoff in the 
GSMF. These indicate that the thermal stellar feedback is 
unable to quench star formation like is observed in massive 
galaxies. 

Comparing SFR with stellar mass, the simulated galax- 
ies lie along a tightly correlated "star forming main se- 
quence." The simulated galaxies match observations by|Ka-| 



jisawa et al. ( 2010 1 at 2 ^ 3, but there is a slight discrepancy 



at z = 2 between simulations and observations. At a given 
stellar mass, the simulated SFRs and correspondingly, the 
sSFRs, are ~ 2 times lower than the observed values at 
9.5 < log(M*/M0) < 10.5. The high sSFRs in low mass 
haloes at z — 2 suggests that there needs to be a signif- 
icant amount of cold gas still present in these galaxies at 
z — 2. Although our model does a better job of delaying the 
star formation at early times than most SAMs and hydro- 
dynamic simulations, after z — 3 the simulated galaxies are 
forming too few stars. 



Dave ( 2008 1 showed that the higher observed SFRs at 



2^2 can be explained by an evolving stellar IMF, which be- 
comes increasingly bottom-light at high redshift. However, 



Marchesini et al. ( 2009 1 showed that when such a bottom 



light IMF was used to model observations, the resulting ob- 
served high-redshift GSMF contained less galaxies, making 
the discrepancy with model GSMFs worse. 

Regarding the evolution of sSFRs at 2 > 3, our simula- 



tion results are consistent with the revised Stark et al. ( 2013 1 



observations for a sample of galaxies with stellar masses cen- 
tred around 5 x lO^M©. The increasing sSFR at z ^ 4 is 
consistent with increasing baryon accretion rates at larger 
redshift translating into larger sSFR in galaxies of fixed stel- 
lar mass. However, our simulated galaxies have lower sSFR 
values than observed ai z = 2. Weinmann et al. (20121 ar- 



gued that the correct sSFR evolution should follow naturally 
from the correct evolution of the GSMF. We see a slight de- 
viation from the observed sSFR relation even though we 
match the GSMF. It must be noted that IWeinmann et alTI 
(20121 performed their analysis at z < 2, while our simu- 
lation has only reached at 2 = 2, where the observational 
estimates are less robust and and might show some internal 
inconsistency among different galaxy properties (e.g. sSFR 
and GSMF). 

There may also be another physical mechanism delay- 



ing star formation. Krumholz et al. (2011 1 and Krumholz & 



Dekel (20121 argue that star formation depends sensitively 



on a metallicity threshold. Until gas reaches this threshold, 
which coincidentally also delays the formation of H2, star 
formation is delayed in low mass galaxies at 2 > 3, which 
leaves sufficient cold gas at 2 = 2 to increase the sSFR of 
these galaxies to the observed values. 

Altogether our results suggest that stellar feedback is 
one of the most important factors regulating star formation 
in Mhaio < 1O^^M0 galaxies. What is most important is 
when the feedback occurs rather than simply the amount of 
feedback energy. Simply increasing and decreasing the feed- 
back energy will only set the normalisation i.e., the total 
stellar mass of present at 2 = 0, but the key is delaying 
star formation in low mass galaxies. When we include stel- 



lar feedback immediately after a star forms until supernovae 
stop exploding after 30 Myr, star formation is significantly 
delayed in low mass galaxies. In this way, we account for the 
downsizing in galaxy populations by delaying the star for- 
mation in low mass galaxies with our stellar feedback model 
and thus reconcile a couple key aspects of a ACDM cosmol- 
ogy with observations. 
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